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ABSTRACT 


In the network diversion problem (NDP), we must find a minimum-weight set of 
edges in a directed graph G = {V,E) whose deletion forces all s-t communication to pass 
through one or more diversion edges in a diversion set E^. We develop and test a 

specialized branch-and-bound algorithm for this NP-complete problem. The algorithm is 
based on partitioning the solution space with respect to edges in certain s-t cuts and yields 
a non-standard, non-binary enumeration tree. The algorithm is coded in Java version 1.4 
and run on a 1.5 MHz Pentium IV computer with 384 megabytes of RAM. An instance 
of NDP on a grid graph with 2502 vertices, 9900 edges and one diversion edge is solved 
in 5.66 seconds; the same problem with 10 diversion edges is solved in only 0.84 
seconds. 
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EXECUTIVE SUMMARY 


In this thesis we eonsider the network diversion problem (NDP) in a direeted 
graph G = (V,E) . The objeetive of NDP is to out links in G, using minimal resouroes, in 
order to foroe all flow from “souroe vertex” s to “sink vertex” t over at least one link in 
speoified set of links. 

NDP arises in the oontext of intelligenoe gathering when trying to foroe an 
opponent to oommunioate over a speoified set of links in a oommunioation network. It 
may also be relevant in warfare, if one wishes to foroe an opponent to travel over 
vulnerable links by destroying links of the opponent’s transportation network. Finally, 
the solution to NDP may also be useful for finding a simple s-t path in a direeted graph 
that oontains a speoified edge (link). 

NDP has been modeled and solved as a binary integer program, but this approaoh 
has not proved very suooessful for problems of even moderate size. 

An s-t out in G is a set of edges whose removal from G disoonneots all direeted 
s-t paths. NDP is stated as the problem of finding a minimum-weight, minimal s-t out 
that oontains a diversion edge. We foous on the problem that has a single diversion edge 
e = (m, v), and use that prooedure as a subroutine to solve the multiple-edge problem. 

In our study, we develop a branoh-and-bound approaoh to NDP that relaxes the 
requirement of minimality in the out identified. This approaoh first identifies a 
minimum-weight s-t out that oontains the diversion edge by “quasi-inoluding” that edge 
in all outs. An edge (m,v) is quasi-inoluded by adding two artifioial, infinite-weight 
edges, one from s to u and one from v to t, so that every finite-weight s-t out in G 
oontains (m,v). If this initial out is a minimal out, NDP is solved. Otherwise, we 
partition the feasible region of NDP using the quasi-inolusion and exclusion technique to 
create an appropriate branching scheme for NDP. An edge is excluded from a cut by 
setting its weight to infinity so that no relevant, finite-weight s-t cut contains that edge. 

We use a non-binary branching scheme where there may be as many branches 
below a node in the enumeration tree as there are edges in a locally-minimum-weight cut. 
The algorithm backtracks whenever a minimal cut is identified or when the weight of the 
locally minimum-weight, non-minimal cut exceeds the weight of the incumbent solution. 
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The new braneh-and-bound algorithm is implemented in the Java 1.4 
programming language, and tested using a 1.5 GHz Pentium IV eomputer with 384 
megabytes of RAM operating under Windows XP. In a 50 by 50 grid network with 2502 
vertiees, and 9900 edges with random weights, an optimal solution to NDP eontaining a 
randomly ehosen link is found in 5.66 seeonds; the same problem with 10 randomly 
ehosen links is solved in 0.84 seeonds. 



I. INTRODUCTION 


The network diversion problem involves eutting links in a network, using the least 
effort possible, to foree the user of the network to eommunieate or travel over eertain 
“exploitable” links. The problem arises in the eontext of intelligenee gathering but may 
also have warfighting applieations. This thesis develops a new, exaet methodology for 
solving this NP-eomplete problem. 

A. PROBLEM DEFINITION AND ITS APPLICATIONS 

We eonsider a directed network (graph) G = {V,E) eomposed of a set of vertices 
V and a set of directed edges E qVxV . Two vertiees are distinguished, a source 
vertex s and a sink vertex t. A speeial diversion set Ej^ cz E is defined; eaeh e g is 
a diversion edge. Eaeh edge e = (m,v)g E has a weight (interdietion eost) w^ > 0 
assoeiated with it; this is the eost or amount of effort required to attaek and destroy the 
edge. We assume that w^ =0 for all e g , and > 0 otherwise. A {directed) path 
in G is a sequenee of edges eonneeting two vertiees, with all edges pointing in the same 
direetion. The network diversion problem (NDP) is; Find a minimum-weight set of 
edges C sueh that all s-t paths in G - C = {V,E\C) must eontain at least one edge of 
Ej^. In other words, if no edges in C ean be used, all flow from s to t must pass through 
one or more diversion edges. 

NDP is stated more eonveniently by defining a minimal s-t cut. This is a minimal 
set of edges C whose removal from G diseonneets all direeted s-t paths. A minimum- 
weight, minimal s-t eut, C ^ E, that eontains at least one edge e& E^ is an optimum 
solution for NDP. 

The intelligenee-gathering applieations of NDP are obvious. However, there may 
also be warlighting applieations in whieh an “interdietor” will interdiet (attaek and 
destroy) links of his opponent’s logisties network and foree that opponent to move 
materiel or troops over eertain links, i.e., diversion edges, that are vulnerable to 
subsequent attaek. 
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A theoretical application of NDP arises in finding a simple (loopless) directed 
path from 5 to t that contains a specified edge ; this problem is shown to be NP- 
complete by Fortune, Hopcroft and Wyllie (1980). A feasible solution to NDP exists, 
with = {e^} (and with =\ for all e & E\ say), if and only if a feasible 
solution to the path-finding problem exists. 

B, BACKGROUND 

NDP has been studied by Curet (2001). He develops an integer-programming 
model (fP) for the problem but finds that fP difficult to solve. He uses heuristics and 
Lagrangian relaxation to solve the problem approximately. The Lagrangian bound could 
be embedded in a branch-and-bound algorithm for solving the problem exactly, but this is 
not implemented. Clearly, the problem is only partially solved; This thesis develops and 
implements a general methodology for solving NDP exactly. 

C. THESIS OUTLINE 

The remainder of this thesis is organized as follows: Chapter If describes our 
general branch-and-bound approach for solving NDP. Chapter Iff explains the algorithm 
in detail. Chapter IV presents computational results and Chapter V summarizes of our 
findings. 
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II. GENERAL METHODOLOGY 


This chapter first reviews basie definitions and notation that will be used 
throughout this thesis. Subsequently, we outline a general branch-and-bound approach 
for solving the Network Diversion Problem. 

A. DEFINITIONS AND NOTATION 

An s-t cut in direeted graph G = {V,E) is a set of edges C such that 
G - C = (V,E \ C) contains no s-t paths. A non-minimal cut is a cut in which one of its 
proper subsets forms another eut. The weight or eapaeity of a cut C is w{C) = Z 

ee C 

The value of a max flow from 5 to t in G equals the eapaeity of a minimum-weight s-t 
cut in G according to the maximum-flow minimum-out theorem (Ford and Fulkerson 
1962, pp. 11-14). 

A rooted (enumeration) tree, T, is a directed oonneoted graph that oontains no 
oyoles, and has a distinguished initial node, oalled the root, whioh does not have any 
inooming edges. There is a simple direeted path from the root node to all other nodes of 
a tree. The nodes on the simple path from a given node i to the root node are ealled 
ancestors of node i. All the nodes that can be reaehed from a node i by going along any 
simple directed path are descendants of node i. The root is said to be on level 0. Any 
other node i in T is at some level / g {1,2,...} which equals the length (number of edges) 
of the path from the root node to node i. 

B, THE NEW BRANCH-AND-BOUND METHODOLOGY 

We propose a completely new braneh-and-bound algorithm for solving NDP, a 
proeedure that is not based on an IP formulation and its linear-programming (LP) 
relaxation. Another unique feature of the methodology is that the assoeiated enumeration 
tree is not a binary tree, but may have as many branches below a node as there are edges 
in an s-t cut. Our description below begins by considering a binary branching scheme, 
however. 
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Consider a typical, binary IP: 


(BIP) 


z* = min cx 

X 

s.t. Ax = b 


X G {0,1}'. 


( 1 ) 

( 2 ) 

(3) 


This is usually solved through its LP relaxation which replaces the integer value 
restrictions (3) with 0 < x < 1: 


(BIPR) z = min cx 

X 

s.t. Ax 

0 < X 



(4) 

= b 

(5) 

< 1. 

(6) 


BIPR is solved forx . If the solution is integer, we have solved BIP; if not, we select 
some fractional Xj and “branch” to solve two problems: 


(BIPRa) z. 


and 


min cx 

X 

(V) 

s.t. Ax = b 

(8) 

0 < X < 1 

(9) 

II 

o 

(10) 


(BIPRb) 


Zj = min cx 

X 

s.t. Ax = b 

0 < X < 1 


( 11 ) 

( 12 ) 

(13) 

(14) 
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An enumeration tree is formed in whieh eaeh problem has the same eonstraints 
and objeetive funetion as BIPR exeept for some additional bounds on certain components 
of X. This branching procedure is repeated recursively, using upper- and lower-bounding 
information, to create a convergent algorithm. The leaves of the tree represent 
subregions of the complete feasible region of BIP. Each subproblem can produce two 
more problems in the enumeration tree when we branch on one of the noninteger 
components of their solution. In order to find an optimal solution to BIP, the 
enumeration tree must be monotonic, i.e., the objective value of each subproblem in the 
tree is no better than that of any of its ancestors. The objective value for the best integer 
solution found so far, called the incumbent solution, is kept as an upper bound on z*. 
Whenever a node is reached whose objective value is no better than the incumbent 
solution, the tree is pruned at that node. The tree can also be pruned at any node if an 
infeasible solution arises from all the restrictions that have been added. In this approach, 
we hope that we quickly find a good integer solution to BIP through one of its restricted 
LP relaxations. 

Now, let us consider an analogous solution procedure for NDP. Assume that 
Ed = ; cases with more than one diversion edge will be discussed later. Let C' 

denote the set of minimal s-t cuts in G, let C" denote the set of non-minimal s-t cuts, and 
let C = C\j C" . We state NDP through this formulation; 


(NDP) 

= min Y, 

ee C 

(15) 


s.t. G C 

(16) 


C G C 

(17) 


This formulation simply says that we are looking for a minimum-weight, minimal s-t cut 
in G that contains the diversion edge . If minimality of the optimal cut C* is not 

enforced, it is possible for C* - to be a cut and thus no s-t communication at all would 
be possible in G - (C* - e^) (Curet 2001). We would have failed to solve the problem 
in this case. 
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We do propose, however, to solve NDP by eonsidering a relaxation of the 
problem in whieh non-minimal euts are allowed; 

(NDPR) z = min ^ 

eeC 

s.t. e C 

C G C' u C" (= C) 

We derive NDPR from NDP by removing the minimality restriction (17), similar to the 
relaxation of BIP to BIPR. We solve NDP using a branch-and-bound procedure based on 
NDPR, analogous to solving BIP through branch-and-bound with BIPR. However, the 
details of branching, bounding, and monotonicity do change. 

We will be solving NDPR through its dual, a max-flow problem. Furthermore, in 
the primal problem, not all cuts of C" are feasible, and the actual set of feasible, non- 
minimal cuts will change as our branch-and-bound algorithm proceeds. However, we are 
always working with valid partitions of subsets of C ", so the following discussion about 
the correctness of the branch-and-bound algorithm is essentially correct. 

The graph G with diversion edge = {u,v) “quasi-included” is denoted 

G @ . The edge is quasi-included in a cut by adding two artificial, infinite-weight 

edges, one from 5 to w and one from v to t. Since any minimum-weight cut must intersect 
at least one of the edges in the path (s,u), (u,v) , and since any minimum-weight 

cut of interest must have finite weight, we are assured that is an element of any 
minimum-weight cut m G ® . 

NDPR can be solved by modifying G into G @ , and by then finding a min- 

weight cut Cq in G © through standard techniques. (That is, find a maximum s-t 
flow in G © using weights as edge capacities, and then identify a minimum cut, in 
©(l^l) time, by finding a cut that is saturated by that maximum flow.) Cq will also be a 
cut of G , and if it is a minimal cut in G , NDP is solved. 


(18) 

(19) 

( 20 ) 
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Unfortunately, the quasi-inclusion method modifies the graph G into G @ so 
that a minimal, min-weight cut in G ® may be a non-minimal cut in G (Balcioglu 
and Wood 2003). Thus, Q = {e^,e 2 ,...,ei^} may not solve our problem. Now, consider 
partitioning the space of feasible solutions based on a particular edge . It will be 
convenient to select from the initial non-minimal cut identified, i.e., e . In 

this case, we create subproblems analogous to BIPRa and BIPRb by alternately including 
or excluding from the initial non-minimal cut: 


(NDPRa) 

2. = min ^ 

ee C 

(25) 


s.t. G C 

(26) 


C G C' u C" 

(27) 


e. e C , 

(28) 


and 


(NDPRb) 

2ft = min ^ 

^ eeC 

(21) 


S.t. G C 

(22) 


C G C' u C" 

(23) 


^ C . 

(24) 


NDPRa and NDPRb represent subproblems in which is quasi-included or 
excluded, respectively. 

NDPRa can be solved using (G ® e^) ® e^. In this case, a new cut C in 
{G ® Cj) ® e^, which is a solution to NDPRa, is found by applying the max-flow min- 
cut theorem. 
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NDPRb can be solved by “exeluding” from all outs of G ® e^. The graph 
G ® with edge exoluded is denoted (G © O . The edge is excluded from 
all outs by setting its weight to infinity. No relevant minimum-weight s-t out oontains e^ 

beoause the solution to NDP must be a out with finite weight. In this ease, a new out C 
in {G @ e^) Q e^, whioh is a solution to NDPRb, is found by applying the max-flow 
min-out theorem. 

NDPRa and NDPRb illustrate a general teohnique, using quasi-inolusion and 
exolusion, to partition the solution spaoe of NDPR. A oonvergent branoh-and-bound 
algorithm oan be oreated through repeated applioations of the partitioning soheme 
analogous to a branoh-and-bound algorithm to solve BIP. Eaoh subsequent node in the 
enumeration tree oorresponds to a restriotion (quasi-inolusion or exolusion of an edge) of 
the problem above it, and thus we maintain monotonioity of the objeotive funotion, whioh 
is required for branoh-and-bound. Monotonioity also derives from the faot that the 
maximum flow in G , whioh equals the oapaoity of the out identified at eaoh node, is 
never deoreasing, beoause quasi-inolusion adds infinite-weight edges and exolusion 
inoreases the oapaoity on an edge. 

The enumeration tree formed by using this partitioning soheme leads to the 
problem of repeatedly generating the same out. NDPRa gives baok the same solution as 
NDPR beoause the edge e^ is already an edge of Q. We use a generalized, non-binary 
partitioning soheme to avoid this. 

Every node in the enumeration tree obtained in the generalized partitioning 
soheme will have as many subproblems as there are edges in the ourrent out. All the 
subproblems below the ourrent non-minimal out Q = {e^,ej^,...,e,^} are derived by 
establishing a partition of the solution spaoe of outs C through the relationship 

C = [C n Uj] u [C n gj n gj] u [C n gj n ^2 ^ 3 ] 

U...U [Cn gj n ^2 G-i ^k\ 

u [C n g; n...n ej. (25) 
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Here, [Cn gj r\ 62 n e^.] denotes the set of all cuts that include through 

, and exclude which is implemented through G © O . 

An optimal solution to NDP must occur in one of the subsets of the partition, 
although it cannot occur in [C n gj n... n e^.], because we have already established that 

any cut in this set must be non-minimal. (Q is non-minimal and any cut in this subset is 
a superset of Q and can thus be ignored.) In other words, including all the edges e\ 
through in a new cut leads to repeated generation of the same cut which we know to 
be non-minimal. 

So, in this instance, our branch-and-bound algorithm will create k branches by 
forcing quasi-inclusion and exclusion of the indicated edges from the initial cut. When a 
minimum-weight, minimal s-t cut contains is found, we can use its weight as a bound 

to trim branches of the tree in the usual manner. A branch is also trimmed at any node if 
an infeasible solution arises from all the edges that have been included in or excluded 
from a current node, i.e., an infinite maximum flow is detected. Note that in solving an 
IP by LP-based branch-and-bound, we allow fractional solutions and “hope for” 
integrality. In our procedure for solving NDP, we allow non-minimality and hope for 
minimality. 
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III. IMPLEMENTATION 


This chapter introduces a simple implementation of the new branch-and-bound 
algorithm for solving NDP, along with a pseudo-code deseription of the algorithm and its 
demonstration on a sample graph. After establishing the correetness of this simple 
implementation, a generalized implementation is stated with a number of refinements. 
Finally, some pathologieal eases and approaehes to their solutions are diseussed. 

A. NEW BRANCH-AND-BOUND ALGORITHM 

This section describes how the new branch-and-bound methodology will be 
implemented to solve NDP. We assume that = {e^} and that we are given a directed 

network. Note that an undirected network can be transformed to a directed network by 
replacing each undirected edge by two directed anti-parallel edges whose weights are 
identical to the weight of the undirected edge (Abuja et al. 1993, pp. 39). 

Algorithm C (Figure 1) to solve NDP is based on the near-min-cut enumeration 
proeedure. Algorithm B, in Baleioglu and Wood (2003). It is stated in simple terms for 
clarity; we discuss a number of praetieal enhaneements later in this ehapter. 


11 



Algorithm C 

INPUT: A directed graph G = {V,E), s,t g V, diversion edge , and edge weights 
w > 0 except that = 0. 

OUTPUT: A min-weight, minimal s-t cut C containing , i.e., a solution to NDP. 

GLOBAL VARIABLES: z ^ oo; C ^ 0 ; 

begin 

<— {e^}; /* set of edges to be included */ 

£“ <— 0; /* set of edges to be excluded */ 

EnumemteC(G , s, t, w, E^ , £“); 

print C; /* If C = 0 , the problem is infeasible */ 

end. 

Procedure EnumerateC{ G ,s,t, w, E^, E~) 

begin 

w' w ; G' G ; 

for (each edge e = {u,v) g E~) w {u,v) <— oo ; 
for (each edge e = (u,v) & E^) begin 

add artificial edge {s,u) to G' and let w'(s,u) <— oo ; 
add artificial edge (v,t) to G' and let w'{v,t) <— oo ; 

endfor; 

/* G' and w' are now interpreted to include artificial edges */ 

/* MaxElowO below finds a min-weight s-t cut Q in G', and its weight z, but subject 

to edges E~ being excluded from the cut, and edges E^ being quasi-included. 

It uses a maximum-flow algorithm to do this. */ 

[z,Co] <— MaxElow{G',s,t,W ); 

if (z > z) return; 

if(Co is a minimal s-t cut in G ) then 

C ^— Cg; z <— z; 

return; 

endif; 

for (each edge bg Cq \ E^) begin 
E <— E u {g} ; 

EnumerateC( G , s,t, w, E^ , E ~); 

if (z > z) return; 

E- ^ E-\ {e}; E^ ^ E^ u {e }; 

endfor; 

return; 

end. 

Figure 1. Algorithm C: A branch-and-bound algorithm to solve the network diversion 
problem when E^ = {e^} . 
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The algorithm starts by assuming no initial solution, i.e., C = 0 and z = oo. The 
diversion edge is added to beeause any solution must contain . Then, the 
algorithm calls the procedure EnumerateC which attempts to find a min-weight s-t cut 
containing . 

The procedure EnumerateC first modifies G into G' by changing the weight of 
each edge being excluded to infinity and by adding two artificial, infinite-weight edges 
for each edge being quasi-included. Then, it finds an initial min-weight s-t cut Cg in G '; 
its weight z is a local lower bound, i.e., a lower bound subject to the restrictions defined 
through E^ and E~. Cg is also an s-t cut in G. If the local lower bound z is greater 

than upper bound z , the algorithm backtracks because no solution better than the 
incumbent can be found given the current restrictions (or any superset of them which 
would occur below in the tree). If a minimal s-t cut whose weight is less than upper 
bound z is found, i.e., a feasible solution to NDP is found, then this cut is kept as an 
incumbent solution, and z is updated. The algorithm also backtracks in this case 
because the weight of any other s-t cut below the current cut in the monotonic 
enumeration tree can be no better than the incumbent solution. In other cases, the 
algorithm calls itself recursively, creating a branch for each edge in Cq\ E^ , i.e., a 

branch for each edge in the current cut yet fixed (forced in). For each branch, it finds a 
new min-weight s-t cut that contains the edges being quasi-included and does not contain 
the edges being excluded according to the generalized partitioning scheme. The 
algorithm terminates when all of the branches are pruned through bounding arguments or 
through infeasibility. The algorithm is clearly finite since the depth of the enumeration 
tree is at most \E\ . 
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Figure 2, Sample Graph. The diversion edge is e^g = (u,v) . All edge weights are 1 

except that w = 0 . 

^ ^10 

To illustrate Algorithm C, consider the enumeration tree in Figure 3, which 
corresponds to solving NDP on the graph of Figure 2. The algorithm first finds an initial 
min-weight s-t cut (node \) m G @ = G @ at the root node at level 

0. This is a non-minimal cut because it contains the minimal cut Thus, the 

algorithm must partition the solution space of cuts based on the edges • 

After partitioning, the algorithm finds a minimal cut at node 2 in 

its first branch marked by ej. It keeps this cut as an incumbent solution, sets z = 4 and 
backtracks (as it always does when it finds a feasible solution). The second branch, 
{ 6 ^, 62 ), yields a non-minimal cut (node 3). The local lower bound is 

z = 3 so the algorithm cannot backtrack. At this node then, it branches to and finds a 

minimal cut (node 4) better than the incumbent. It keeps this cut as the incumbent 
solution and backtracks while setting z = 3. Nodes 5 and 6 are pruned by bounding, 
without exploration, because the corresponding local lower bounds cannot be better than 
z , i.e., no recursive call is made for nodes 5 and 6 because of the third return statement. 
Finally, the algorithm terminates with the optimal solution C* = . 
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B, CORRECTNESS OF THE ALGORITHM 

To state that Algorithm C correctly solves the Network Diversion Problem, it 
suffices to show that (a) the algorithm would enumerate all feasible cuts containing the 
specified diversion edge if backtracking occurs only by infeasibility, and (b) the 
enumeration tree maintains monotonicity of the objective function. 

Balcioglu and Wood (2003) give a proof that all s-t cuts in a graph can be found 
by using this partitioning scheme with an initial min-weight s-t cut. In our methodology, 
we begin to partition by quasi-including the diversion at the root node of the enumeration 
tree; but this is the same technique that Balcioglu and Wood (2003) use in their 
enumeration algorithm, to force inclusion of edges. Hence, Algorithm C with no “value- 
based backtracking” must enumerate all minimal cuts containing . 

When the solution space is partitioned recursively, the weight of an edge in G 
has been increased to infinity or two infinite-weight edges have been added to G at each 
node of the enumeration tree. The weight of a cut below a current node never decreases 
because this process never reduces the maximum flow in G. Therefore, we conclude that 
the enumeration tree is mono tonic. 
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C. REFINEMENTS TO ALGORITHM C 

Algorithm C can be refined in a number of ways for improved effieiency and 
greater generality: 

1. Algorithm C, as stated above, repeats mueh work for the sake of elarity, but 
the aetual implementation avoids this for effieiency. For instance, the implementation 
does not modify the original graph from serateh at eaeh node of the enumeration tree (by 
using and E~ ). Rather, onee an edge is quasi-ineluded or excluded from all cuts at 
some node in the enumeration tree, it is automatieally quasi-ineluded or exeluded from all 
deseendant-node problems, with no additional work required. 

2. A max flow in a parent-node graph is always feasible for all its ehild-node 
graphs. This makes identifieation of a min-eut in the ehild nodes signifieantly faster, 
through a small number of flow augmentations, because we do not have to identify a max 
flow and thus a min cut from scratch. 

3. Most branch-and-bound algorithms can be improved by adding a heuristic that 
finds a feasible initial solution so that the algorithm begins with z < oo. Nodes of the 
enumeration tree ean often be trimmed quickly this way, rather than having to wait for 
the algorithm to stumble upon a feasible solution and eorresponding finite upper bound. 
The problem of finding a minimal cut containing a specific edge is an NP-eomplete 
problem (Caryle and Wood 2002). Thus, the following heuristie is not guaranteed to 
work (Curet 2001); in praetiee it usually does, however. Let = (m,v) : 

a. Using breadth-first seareh (BFS), find a shortest path (minimum number 
of edges) from s to u that traverses neither v nor t. Let E^ be the set of 
edges in that path. 

b. Using BFS, attempt to find a shortest path from v to t that does not use any 
vertex on the first shortest path. If sueh a path is found, let its edges be 
denoted E^. 

e. If two paths ean be found, a min-weight eut in G O {E^ u E’j) will be a 
minimal eut eontaining (Curet 2001). 
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We implement the heuristic as a subroutine TwoPathsHeuristic(G, s, u, v, t, w) 
which returns a finite upper bound z if successful, and returns 0 ° otherwise. Figure 4 
gives an example. 



Figure 4. Illustrating the TwoPathsHeuristicO to find a feasible solution for NDP. A 
shortest s-u path is first identified as = {e^,e^} . Then, a shortest v-t path, which does 

not use any vertex on the first shortest path, is identified as = {e^ 2 ,eg}. Finally, a 
minimal cut C' = {e 2 ,e 4 ,e^,e^,eig,eii} is computed in G O (E^ vj E 2 ), and w{C') = 5 is 
kept as an initial upper bound for NDP when = {e^g}. 

4. We may wish to solve NDP when = {ej,...,e^}, i.e., when |£'^| > 1. This 
is easily accomplished by running Algorithm C with = e^, then with = 62 , etc., and 

taking the best solution from among the p solutions. However, it is clear that we can 
improve that procedure by using bounding information and partial solutions from one 
subproblem to the next. Furthermore, if \Ej 2 | is large, it may be worthwhile to first find a 

min-weight cut C' in G and simply check to see if some diversion edge is contained in 
that cut. (Recall that = 0 for all e g , and all other edges have positive weights. 

Thus, if \Ej 2 I is large, it seems intuitively likely that one of the edges in will turn up 
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in a min-weight cut.) If it does, NDP has been solved by just a single min-cut (max-flow) 
calculation. If not, the flow just found is still a useful point from which to begin the 
algorithm. 

The improved and generalized implementation of Algorithm C to solve NDP when 
|£'^| > 1 is shown in Figure 5. 


Algorithm C 

INPUT: A directed graph G = {V,E), s,t € V, diversion edges and edge 
weights w > 0 except that w = 0 for all e g . 

OUTPUT: A min-weight, minimal s-t cut containing at least one e i.e., a 
solution to NDP. 

GLOBAL VARIABLES: z ^ oo; C ^ 0 ; 

begin 

[w(Co),Cq] <— MaxElow{G,s,t,yi) 
if(Co n ^ 0) then print Cg and halt; 
for (each edge e g ) begin 

^ {^d) i 

u <— tail vertex of ; v <— head vertex of ; 

[z,Cg] <— TwoPathsHeuristic{G,s,u,v,t,'w); 

if (z < z ) then 

z 4— z C ^— Cg; 

endif; 

endfor; 

for (each edge e g ) begin 

E^ <— ; /* set of edges to be included */ 

E” <— 0 ; /* set of edges to be excluded */ 

EnumerateC( G , s,t, w, E '^, E~ ); 

endfor; 

print C; /* If C = 0 , the problem is infeasible */ 

end. 

Procedure EnumerateC{ G ,s,t, w, E^, E~) 

begin 

Same as in Figure 1 ; 

end. 

Figure 5. An improved and generalized version of Algorithm C for \Ej^\ > 1. 
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D, PATHOLOGICAL CASES 


We call an edge that is not contained in any simple s-t path in G an extraneous 
edge. A pathological case arises when all diversion edges are extraneous edges for NDP; 
In this case, a solution to NDP does not exist but Algorithm C “solves” the problem 

correctly and returns C = 0 . However, the algorithm may perform a huge number of 
iterations in such cases. Pathological cases cannot occur, however, if the extraneous 
edges in a network are identified and marked as “ineligible” to be diversion edges. 

Before examining pathological cases, we need to define some concepts associated 
with undirected graphs. Any vertex in a connected undirected graph whose removal 
results in a disconnected graph is called a separation vertex. A graph containing a 
separation vertex is called separable, and a graph that does not contain such a vertex is 
called nonseparable. A subgraph G' of G that contains no separation vertices is a 
nonseparable component of the original graph. 

Even (1979, pp. 57-62) gives an algorithm, based on depth-first search (DFS), for 
finding all the nonseparable components and the separation vertices of a given graph in 
0(|E'|) time. Consider the separable graph G in Figure 6 to find its nonseparable 
components. 



Figure 6, A separable graph G with separation vertices b, e, and i. 

First, we modify G into G -I- (5,t) by adding an artificial edge from 5 to t in 
order to keep s and t in the same nonseparable component. See Figure 7. The 
nonseparable components (Figure 8) are then found with Even’s algorithm. 
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Figure 7. The graph of Figure 6 modified to G + {s,t). 



Figure 8. The nonseparable components of the graph G + {s,t) from Figure 7. 

The extraneous edges in the graph can now be determined easily: All edges of 
any separable component of G + {s,t) , except for the component containing both s and t, 

are extraneous. For instance, the nonseparable components G^ and Gj of G + ( 5 ,t) do 
not contain both 5 and t. Therefore, the edges (e,/), {i,g) ,{g,h), and {h,i) from Gj 
and G 3 are extraneous, and are not allowed to be diversion edges. Thus, we are 
effectively operating on the reduced graph G" of Figure 9. 
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Figure 9, The reduced graph G" after nonseparable decomposition. The artificial edge 
{s,t) is removed because it is not useful any more. 

The nonseparable decomposition process applied to undirected graphs is 
guaranteed to identify all extraneous edges, and this process is also valid for finding 
extraneous edges of a directed graph if it is treated as being undirected. But the 
nonseparable decomposition process may not identify all extraneous edges in a directed 
graph. In fact, there is no known polynomial-time algorithm to detect all extraneous 
edges in a directed graph (Fortune, Flopcroft and Wyllie 1980). This is the same NP- 
complete problem as finding a simple directed path from 5 to t that contains a specified 
edge. Therefore, extraneous edges in large directed graphs may cause Algorithm C some 
difficulties. 

Our algorithm does not include a nonseparable decomposition procedure because 
the networks that we test do not contain any nonseparable components. But it may be 
worthwhile to implement this process because many real-world networks may contain 
extraneous edges. 

If an undirected network is being analyzed, we do not know in which direction an 
interdiction edge = (m,v) should appear in the optimal solution. That is, we do not 

know if all s-t paths in G - C* -I- (w, v) traverse from w to v or from v to u. We can solve 
the problem, in theory, by quasi-including the directed edge = {u,v) and running the 
algorithm, then quasi-including the directed edge = (v,m) and running the algorithm 
again, and then selecting the best of the two solutions as the optimal solution to NDP. 
However, it is possible that the “optimal solution” for one of the orientations is C* = 0 . 
In this case, the algorithm may fail to find an optimal solution in a reasonable amount of 
time for large networks. There is a simple way to deal with this issue. 
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Recall that an undirected edge is represented by two anti-parallel directed edges in 
Algorithm C. Figure 10 shows the representation of G" from Figure 9, which would be 
used in Algorithm C. 



Figure 10. The directed representation of G" (Figure 9) used in Algorithm C. 

Let us refer to the undirected edge (/, k) in G' as {/, k]. The flow from 5 to t in 
the directed representation of G" travels over either {i,k) or {k,i)- The removal of 
either (/, k) or {k, i) represents that the edge {/, k) in G" is in a cut. If the edge (/, k) is 
chosen to represent the diversion edge, then an optimum solution exists for NDP. On the 
other hand, if the edge {k,i) is chosen, then no feasible solution is obtained to NDP: 

Algorithm C might fail to return the answer C* = 0 in a reasonable amount of time for 
larger networks with similar structure. 

To avoid this, we replace the edge {i,k} with a directed Wheatstone bridge as 
shown in Figure 11. The edge (v^Vj) that takes place of the edge {i,k} with identical 
weight is assigned as the diversion edge, i.e., two artificial, infinite-weight edges ( 5 ,Vj) 
and (Vj,^) are added to the directed representation of G" . In this case, all finite-weight 
cuts in the new graph must contain the edge (Vj,V 2 ). Because of that, we conclude the 
edge {/, k) always occurs in a cut oriented in whichever direction is appropriate. 
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Figure 11. Representation of the diversion edge {i,k} in G" using a Wheatstone 
bridge. 
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IV. COMPUTATIONAL RESULTS 


This chapter provides computational results for Algorithm C on a set of structured 
test problems. The algorithm is written using the Java 1.4 programming language. All 
tests are performed on a personal computer with a 1.5 GHz Pentium IV proeessor, 384 
MB of RAM and running under the Windows XP operating system. 

A. TEST NETWORKS 

Topology is the physieal eonfiguration of conneetions eomprising a network, i.e., 
it is the shape of the network. We use two eommon topologies in eommunieations 
networks to test Algorithm C. These topologies ensure high reliability beeause of their 
high. 

1. Grid Networks 

A grid topology consists of a network where the vertices may be viewed as being 
points on the LxH integer lattice in the non-negative quadrant. Each vertex is eonneeted 
to other vertices, which are positioned in the lattice directly to its north, south, east, and 
west. There are two additional vertices, a souree and a sink vertex, which are positioned 
at the west and east side of the network, respeetively. All westernmost vertiees are 
eonneeted to the souree vertex, and all easternmost vertiees are eonneeted to the sink 
vertex. The number of vertiees in a grid network is |E| = HxL-I-2, and the number of 

edges is |E'| = 4 xHxL-2xL. Edges incident to s and t have infinite weights and eannot 

be interdieted and cannot be diversion edges. Eigure 12 shows an example of a small 
grid network. 
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Figure 12. An unweighted, undireeted grid network with H=3 and L=4. The 
undirected edges between two vertices represent two anti-parallel directed edges, and 
their weights are 1. 


A grid topology, or something close to it, is often used in Wide Area Networks 
(WANs) where reliability is important (Harris 1999). 

1. Star-Mesh Networks 

A star-mesh network (Miller 2001) is created if H “rays” of L vertices, each 
connected by an additional center vertex, are cross-connected as well as radially 
connected. The source vertex s is the center vertex. The sink vertex t is chosen from the 
vertices which are positioned at the outer most ring of the network to increase the 
distance between s and t as much as possible. The number of vertices in a star-mesh 
network is |F| = HxL-I-1, and the number of edges is |£'| = 4xHxL . Figure 13 shows an 

example of a star-mesh network. Star-mesh networks are often used in local area 
networks (LANs) (Sharma 1990, pp. 9-12) where the central vertex controls all the 
vertices in the network. 



Figure 13. An undirected star-mesh network with H = 8 and L = 2 . 
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B, RESULTS 

We use a grid network generator (GGFGEN) written in Java by Balcioglu and 
Wood (2003) to generate grid networks. We have coded a star-mesh network generator 
(STARGEN) in Java to generate star-mesh networks. The input parameters of the 
generators, H and E, determine the size of the generated network. 

We test Algorithm C on unweighted and weighted networks separately. A 
network is “unweighted” if all interdictable edges e have = 1. In our weighted 
networks, for each interdictable edge is randomly drawn from the discrete uniform 
distribution on [l,5]. Eor both weighted and unweighted grid networks, the edges 
incident to s and t have infinite weights. 

The sizes of networks generated are systematically increased to observe the 
performance of the algorithm. Each network is also tested as \Ej^\ increases; the set Ej^ 

is chosen randomly from E. 

1. Unweighted Networks 

Table 1 presents problem and solution statistics for Algorithm C applied to 
solving NDP in unweighted, undirected grid networks. Eor all networks in the table, the 
algorithm gives an optimal solution in less than one second; The TwoPathHeuristicQ 
always finds the optimal solution, and this solution is proven optimal at the root node of 
the enumeration tree. Indeed, with some thought about the structure of these problems, 
we can see that they can be solved by inspection. 

Table 2 presents processing times on unweighted, undirected star-mesh networks. 
In this network type, processing times, upper bounds, and optimal solutions change 
suddenly depending on the diversion edges and sink vertex which are chosen randomly. 
Increased numbers of non-minimal cuts in the enumeration tree imply longer processing 
times. In fact, some of these problems, and the weighted problems for star-mesh 
networks, cannot be solved to optimality in a reasonable amount of time. In these cases, 
we have identified s-optimal solutions, for rather large values of s > 0, by replacing 
both instances of “if (z > z) then” with “if (z (1 -I- s) > z) then.” 
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Table 1. Problem statistics and computational results for Algorithm C on unweighted, 
undirected grid networks. GRIDNET HxL denotes a grid network with H rows of L 
vertices. zUB is the upper bound obtained from the TwoPathsHeuristicQ procedure. zLB 
shows the weight of the cut at the root node of the enumeration tree. The column s “Number 
Non-minimal Cuts” and “Number Minimal Cuts” give the number of non-minimal cuts and 
minimal cuts in the enumeration tree, respectively. Here, the algorithm does not need to 
perform any enumeration, and therefore these entries are all zero. The total number of nodes 
in the enumeration tree is the sum of these numbers plus one. 
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Network 


1^1 

\Ed\ 

zUB 

zLB 

w(C*) 

Number 

Non-minimal 

Cuts 

Number 

Miuimal 

Cuts 

Total Time 
(sec.) 




1 

6 

6 

6 

0 

0 

0.00 




2 

8 

6 

8 

6 

0 

0.02 

STARNET 



3 

11 

6 

11 

234 

0 

0.08 


101 

400 








10x10 



4 

7 

6 

7 

4 

0 

0.02 




5 

10 

6 

10 

132 

0 

0.05 




10 

4 

4 

4 

0 

0 

0.02 




1 

2 

2 

2 

0 

0 

0.02 




2 

10 

6 

10 

74 

0 

0.17 

STARNET 



3 

14 

6 

14 

3963 

0 

6.45 


401 

1600 








20x20 



4 

17 

6 

17 

110608 

0 

197.30 




5 

13 

6 

13 

3456 

0 

6.16 




10 

9 

6 

9 

107 

0 

0.25 




1 

29 

6 

29 (0.75)* 

27449 

0 

170.39 




2 

11 

6 

11 

163 

0 

0.81 

STARNET 

901 

3600 

3 

27 

6 

27 (0.50)* 

210612 

0 

1068.39 

30x30 



4 

29 

6 

16 (0.75)* 

83083 

1 

455.78 




5 

16 

6 

16 

56305 

0 

298.16 




10 

23 

6 

23 (0.25)* 

102270 

0 

551.97 




1 

39 

6 

39 (1.25)* 

70205 

0 

710.28 




2 

10 

6 

10 

74 

0 

0.87 

STARNET 



3 

21 

6 

21 (0.25)* 

113284 

0 

1173.43 


1601 

6400 








40x40 



4 

30 

6 

30 (1.00) 

15569 

0 

177.23 




5 

21 

5 

21 (0.50)* 

195531 

0 

2050.08 




10 

23 

6 

23 (0.75)* 

14551 

0 

154.69 




1 

26 

6 

26 (0.50)* 

3960 

0 

74.22 




2 

39 

6 

39 (1.50)* 

20212 

0 

344.90 

STARNET 



3 

21 

6 

21 (0.50)* 

4960 

0 

84.83 


2501 

10000 








50x50 



4 

28 

6 

28(1.00) 

5330 

0 

91.27 




5 

49 

6 

49 (2.50)* 

7190 

0 

127.28 




10 

3 

3 

3 

0 

0 

0.84 


* This problem could not be solved to optimality in a reasonable amount of time, so an s-optimal solution is given. 
The value of e is given in parentheses. 


Table 2, Problem statistics and computational results for Algorithm C on unweighted, undirected 
star-mesh networks. STARNET HxL denotes a star-mesh network with H rays of L vertices 
each. The difference between zUB and zLB is much larger for the more difficult problems. 
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2. Weighted Networks 

Tables 3 and 4 present problem and solution statisties on weighted, undireeted 
grid and star-mesh networks, respeetively. We use the same diversion sets on 
unweighted and weighted networks. Integer edge weights are generated randomly on 
[ 1 ^ 5 ]. 









Number 

Number 

Total Time 
(sec.) 

Network 

I'd 

1^1 

I^dI 

zUB 

zLB 

w(C*) 

Non-minimal 

Minimal 








Cuts 

Cuts 




1 

38 

28 

35 

146 

1 

0.06 




2 

22 

22 

22 

0 

0 

0.02 

GRIDNET 

102 

380 

3 

20 

20 

20 

0 

0 

0.02 

10x10 



4 

24 

24 

24 

0 

0 

0.02 




5 

26 

26 

26 

0 

0 

0.03 




10 

19 

19 

19 

0 

0 

0.03 




1 

48 

42 

48 

1913 

0 

4.38 




2 

48 

48 

48 

0 

0 

0.03 

GRIDNET 

402 

1560 

3 

48 

43 

48 

635 

0 

1.48 

20x20 

4 

37 

37 

37 

0 

0 

0.05 




5 

37 

37 

37 

0 

0 

0.06 




10 

39 

39 

39 

0 

0 

0.06 




1 

82 

80 

82 

62 

0 

1.02 




2 

79 

77 

79 

48 

0 

0.83 

GRIDNET 

902 

3540 

3 

73 

73 

73 

0 

0 

0.08 

30x30 

4 

75 

75 

75 

0 

0 

0.09 




5 

76 

76 

76 

0 

0 

0.09 




10 

73 

73 

73 

0 

0 

0.11 




1 

84 

84 

84 

0 

0 

0.13 




2 

89 

89 

89 

0 

0 

0.14 

GRIDNET 

1602 

6320 

3 

88 

93 

88 

0 

0 

0.16 

40x40 

4 

90 

89 

90 

8 

0 

0.81 




5 

87 

87 

87 

0 

0 

0.17 




10 

85 

85 

85 

0 

0 

0.20 




1 

123 

121 

123 

50 

0 

5.66 




2 

122 

121 

122 

28 

0 

3.05 

GRIDNET 

2502 

9900 

3 

112 

112 

112 

0 

0 

0.25 

50x50 



4 

112 

112 

112 

0 

0 

0.25 




5 

111 

111 

111 

0 

0 

0.28 




10 

119 

117 

118 

3 

1 

0.84 


Table 3. Problem statisties and eomputational results for Algorithm C on weighted, 
undireeted grid networks. The problems that have some non-minimal cuts and no minimal 
cut show that the feasible solution (zUB) found in the TwoPathsHeuristicO is optimal. The 
problems with some non-minimal and minimal cuts show that the optimal solution is found 
in the enumeration tree. 
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Network 


1^1 

\Ed\ 

zUB 

zLB 

w(C*) 

Number 

Non-minimal 

Cuts 

Number 

Miuimal 

Cuts 

Total Time 
(sec.) 




1 

23 

19 

23 

4 

0 

0.02 




2 

35 

21 

32 

131 

1 

0.06 

STARNET 



3 

35 

15 

21 

29 

3 

0.03 


101 

400 








10x10 



4 

12 

12 

12 

0 

0 

0.02 




5 

25 

16 

19 

21 

2 

0.03 




10 

9 

9 

9 

0 

0 

0.02 




1 

13 

13 

13 

0 

0 

0.02 




2 

35 

20 

35 

926 

0 

2.19 

STARNET 



3 

29 

14 

29 

405 

0 

0.91 


401 

1600 








20x20 



4 

41 

15 

36 

21139 

3 

47.09 




5 

24 

17 

24 

49 

0 

0.16 




10 

28 

14 

25 

675 

2 

1.42 




1 

69 

21 

69 (0.50)* 

19068 

0 

137.75 




2 

21 

18 

21 

6 

0 

0.17 

STARNET 

901 

3600 

3 

66 

24 

66 (0.50)* 

19001 

0 

153.12 

30x30 



4 

61 

18 

40 (0.50)* 

27971 

1 

203.89 




5 

38 

16 

38 

8771 

0 

40.34 




10 

56 

16 

56 (0.50)* 

18711 

0 

135.36 




1 

97 

16 

97 (1.25)* 

19978 

0 

315.90 




2 

21 

16 

21 

9 

0 

0.33 

STARNET 



3 

49 

13 

49 (0.25)* 

62458 

0 

925.06 


1601 

6400 








40x40 



4 

87 

11 

87 (1.25) 

23710 

0 

376.97 




5 

62 

15 

62 (0.50)* 

53018 

0 

829.11 




10 

50 

10 

50 (0.50)* 

20310 

0 

315.41 




1 

115 

21 

115 (1.50)* 

8774 

0 

203.16 




2 

118 

17 

118 (2.00)* 

5090 

0 

127.12 

STARNET 



3 

54 

21 

54 (0.25)* 

26486 

0 

673.65 


2501 

10000 








50x50 



4 

84 

16 

84(1.00) 

21363 

0 

513.34 




5 

118 

18 

118(2.00)* 

13662 

0 

365.40 




10 

57 

15 

57 (0.50)* 

11930 

0 

322.52 


* This problem could not be solved to optimality in a reasonable amount of time, so an s-optimal solution is given. 
The value of s is given in parentheses. 


Table 4, Problem statistics and computational results for Algorithm C on weighted, undirected 
star-mesh networks. 
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c. 


INTERPRETATION 


In the grid networks, the differenee between the weight of the feasible solution 
found by TwoPathHeuristicQ (upper bound) and the weight of the solution found at the 
root node of the enumeration tree (the minimum-weight eut in G © ) is small so that 

the algorithm does not need to enumerate many euts. In faet, G @ is likely to yield an 
optimal solution in the grid networks because s and t are highly connected. 

In the star-mesh networks, the typical initial minimum-weight cut from G @ is 

the non-minimal cut that cuts edges incident to t, plus the diversion edge. But the optimal 
solution will usually have many more edges in it. There is a large gap between the upper 
bound and the initial cut. Because of this, the algorithm must enumerate a large number 
of cuts. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The network diversion problem (NDP) arises when a set of links in a network 
must be removed using the least effort in order to foree the user of the network to 
eommunieate or travel over a speeified set of links between two designated vertiees. This 
NP-eomplete problem may have intelligenee-gathering and warfighting applieations. In 
this thesis, we have developed and implemented a general braneh-and-bound 
methodology to solve NDP. 

NDP can be modeled as a binary integer problem (IP), but the literature indicates 
that these IPs are difficult to solve. Therefore, we introduce a new methodology, a 
general braneh-and-bound approach, to solve NDP. The methodology is not based on 
solving linear-programming relaxations. Rather, the minimality of the cuts identified is 
relaxed. In this methodology, the algorithm first finds a minimum-weight s-t cut 
containing the diversion edge which is quasi-included. If this initial cut is non-minimal, 
then a monotonic enumeration tree is formed by forcing quasi-inclusion and exclusion of 
edges from the current cut. An edge (m,v) is quasi-included in a cut by adding two 
artificial, infinite-weight edges, one from s io u and one from v to t. An edge is excluded 
from a cut by setting its weight to infinity. The associated enumeration tree is non-binary 
since it may have as many branches below a node as there are edges in an s-t cut. 

The algorithm has been improved by incorporating certain internal efficiencies, 
and by adding a heuristic to find a good initial feasible solution. When there is more than 
one diversion edge, the solution to the problem can be found by solving NDP for each 
diversion edge separately, and taking the best of those solutions. However, by using 
bounding information and partial solutions from one subproblem to the next, faster 
solutions are obtained. In fact, it is often easier to find a solution when the number of 
diversion edges is large. 

We have tested our algorithm on two common communication-network 
topologies, grid and star-mesh topologies; with diversion edges chosen randomly. We 
increased the size of networks and the size of the diversion-edge set systematically to 
observe the performance of the algorithm. Grid networks were tested with up to 2502 
vertices and 9900 edges (a 50 by 50 grid network). An optimal solution to NDP for all 
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sizes of grid networks is obtained in seconds. For example, in an unweighted grid 
network with 2502 vertices and 9900 edges (a) the algorithm returns an optimal solution 
in 0.125 seconds when \Ejj\ = 1, and (b) it takes only 0.172 seconds to get an optimal 

solution when | = 5 . 

We have also tested star-mesh networks with up to 2501 vertices and 10000 edges 
(a 50 by 50 star-mesh network). An optimal solution to NDP is obtained in a reasonable 
amount of time for networks as large as 20 by 20. Only some of the larger networks can 
be solved optimally. 

Although the Algorithm C works efficiently for certain test networks, its overall 
run time and efficiency can be improved in different ways. For example, if the number of 
non-minimal cuts in an enumeration tree of a NDP can be decreased, the algorithm will 
be more efficient, and run faster. To accomplish this, the quasi-inclusion technique that 
creates non-minimal cuts should be investigated for improvement or replacement. 

Another practical improvement for computation times on some large networks 
can no doubt be obtained by solving the initial maximum-flow problem using a more 
efficient algorithm, e.g., binary blocking flow algorithm (Goldberg and Rao 1998). 

For problems like the star-mesh networks, we may wish to partition based on 
some set of edges other than those in the current minimum-weight cut. For instance, 
branching on edges in a path from the diversion edge to 5 or t may prove more effective. 
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